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Abstract 


Surface roughness spectra of nine radar backscatter units in the 
Askja Caldera region of Iceland were predicted from computer-enhanced 
like- and cross-polarized radar images. A field survey of the caldera 
was then undert tken to check the accuracy of the preliminary anaiysis. 
There was good agreement between predicted surface roughness of back- 
scatter units and surface roughness observed in the field. In some 
cases, variations in surface roughness could be correlated with pre- 
viously mapped geologic units. 
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SECTION I 
INTRODUCTION 


Until recently, side-looking airborne radar images have been used 
in geological investigations in a similar way to aerial photography to 
extract morphologic and structural information. Although viewing geom- 
etry is different from that of conventional aerial photography, many of 
the same photointerpretation techniques can be used in the interpretation 
of radar data. Radar images clearly show topography, linear features, 
fault scarps, and drainage patterns. A second class of information can 
also be obtained from radar data since the brightness in each resolution 
cell of the Image is related to the surface roughness and dielectric con- 
stant of the imaged terrain. 

Recent studies in Death Valley by Schaber et al. (197f>) and Daily 
et al. (1978) have demonstrated correlations between radar reflectivity 
and the lithology of the surface in alluvial deposits. The purpose of 
this study was to verify current interpretation techniques such as those 
used in the Death Valley studies and to investigate the extension of 
these techniques to nonalluvial terrain. This study involved a test of 
radar image interpretation which included an analysis of computer- 
enhanced radar images to determine surface roughness and a subsequent 
field check to evaluate the accuracy of predictions. The Askja Caldera 
region in north-central Iceland was chosen as a test site because of 
the diversity of volcanic surface textures and a paucity of vegetation 
in the area. 

A brief review of the geological setting of Askja is presented in 
Section II. A general discussion of radar image properties is included 
in Section III. Computer processing techniques used to enhance the 
radar data are discussed in Section IV. The correlations between the 
radar imagery and field observations are described in Section V. 
Conclusions are presented in Section VI. The theoretical basis for 
radar backseat ter studies is discussed in the Appendix. 
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GEOLOGIC S KITING 


Miocene plateau basalts, dated at lb million years BP, are the 
oldest rocks in Iceland (.Beloussov and Milanovsky, 1977). Between 
Middle and Late Miocene times, approximately 9 million years ago, volcanic 
activity occurr* d in northeastern and west central Iceland. Subsidence 
ot the middle part ot Iceland that r suited in a slight downwarping of 
the plateau basalts probably also <; urred about this lime. Gradual 
subsidence of the crust changed to tracturing and block faulting during 
a third stage of volcanic activity, approximately 3 million years ago. 
Volcanism was concentrated in southern Iceland in a band approximately 
250 km wide. In northern Iceland, volcanic activity occurred in two 
areas separated by a stable massif. Basalts and palagonit izcd basaltic 
glass were deposited during this stage. 


The central grabens of the modern rift system of Iceland were 
formed approximately 700,000 years ago. The median zone dissects 
Tertiary Plateau basalts from southwest to northeast. In the south, the 
zone i separated into two segments by an Upper Tertiary-Quaternary 
anticline dig. 1). Approximately one-third of the lava produced on 
earth during the past 500 yearu was erupted in the median zone of Iceland 
(Thorar insson , 1965). 
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Fig. 1. Geologic map of Iceland (from Williams et al., 1973) 


Askja Caldera, an active volcanic center, is located within the 
median zone in north-central Iceland. The caldera lies within the 
Dyngjufjoil massif at approximately 65 03 'N, 16 55 'W. The massif con- 
sists of late Pleistocene and Holocene palagonite tuffs and breccias and 
basaltic Lavas (Thorar insson and Sigvaldason, 1962). Early eruptions of 
Askja, centered along fissures near the present-day caldera rim, 
accompanied the postglacial uplift of Dyngjufjoil. Some of the early 
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lavas formed portions ol the Otbrunl lava t leld to the north and the 
Vikursandur lava : leld to the east of Askji. Other lavas flowed south 
and termed a plain between Askja and the nyngjujttkul 1 Icecap. These 
lavas also formed the t loor ot tin precaldera cauldron (Bemmelen ar.d 
Rut ten, 1955). The Askia Caldera w a> formed during tne final stages of 
deformation of Dyngjuf J811 . Basaltic lava poured out from the marginal 
fractures at the western and southwestern »ides ot Askja and covered 
most ot the caldera t loor after caldera subsidence. 

Askja has erupted four times since 1870. In 1875, a Large pyro- 
clastic eruption tormed the crater V £ 1 1 in the northeastern corner of the 
caldera. After the 1875 eruption, caldera subsidence created Oskjuvatn, 
which is 11 km _> In area (Thorar insson and Sigvaldason. 1962). Basaltic 
eruptions occurred twice between 1921 and 1929, both within and outside 
of the caldera. First, four small flows were produced between 1921 and 
1921: Bat shraun (1921), Mvvetningahraun (1922). and two small flows In 

the southeast corner of Askja (1921); then, between 1924 and 1929, the 
large flow located south of Askja poured out from a 6-km-long fissure on 
the southern flank of the Dyngjuf |flll massif. In 1961, eruptions north 
of the caldera produced a basaltic lava that flowed eastward through a 
break In the caldera called ttsk)uop (Thorar insson and Sigvaldason, 1962; 
see Fig. 2). 
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Fig. 2. Geologic map showing historic fl'* in Askja region 
(modified from Bemmelen and Kutten, 1955) 
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SECTION III 
RADAR IMAGES 


Hu* Askj.i Cal di • » region wan chosen for a radar tost site 
because of the diversity of surface textures and the paucity of vegeta- 
tlon. The radar Images used In this study were obtained In July 197b 
with a 25-cm-wavelength, side-looking airborne radar (see Table 1 tor 
radar system specification). In this study, horizontally polarized 
radiation was transmitted and horizontally and vertically polarized 
backscutter were recorded. These modes are referred to as HH and HV, 
respectively. The radar equipment was flown aboard a NASA Convalr 990 
Jet aircraft operated by the Airborne Sciences Office of Ames Research 
Center. The images were obtained at an altitude of 10,000 m. The look 
angle across the images varies from 0 to 45 deg. The flight line is 
shown in Fig. 3 

Table 1. 


Ope rat ing t requency 
Ope rat in, , wavelength 


1. radar system specifications 


1.2 Ghz 
25 cm 




\ 

I 


Polarization Horizontal transmit, 

horizontal receive; 
Horizontal transmit, 
vertical receive 


Azimuth resolution 25 m 

Range resolution, m 15/sin 0 , where is the 

Incidence angle which varies 
from 6 ■ 0 to 0 ■ 45 deg 



Fig. 3. Location map of Askja Caldera showing swath covered by radar 

Synthetic Aperture Radar uses coherent microwaves to illuminate 
the surface below and to the side of the aircraft. The radar receiver 
detects the echo returned from the surface. The returned signal is 
mixed with a reference signal, and the resulting inferterence signal is 
recorded optically on film. The resulting signal film is developed on 



I . rottl I .nul 1 1 1 urn inated with coherent light in m I'pt ii il correlator* 
In thin process, a visible-light replica ol the radar waves received by 
the antenna Is formed that recreates an image of tin. original terrain 
(see e.g., Jensen et al., 1 V 7 7 . for a detailed discussion). Radar Image 
t l 1m can be printed as photographs oi digitized for computer processing. 
Kxamples of HH ( 1 ike-polar Ized; and HV (cross-polarized) Images are shown 
In Figs. 4 and a. 
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Like-polarized (HH) radar image of Askja Caldera 
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Fig. 5. Cross-polarized (HV) radar image of Askja Caldera 

Hie bright line at the top of the Images Is the return from directly 
'•clow the aircraft (nadir), where the look angle ■ 0 deg. Geometric dis- 
tort ion^ in the images Include compression of the i age near nadir, radar 
foreshortening, and shadowing. Radar shadows occur on the side of a 
terrain feature most distant from the transmitter because the radiated 
energy does not pass through objects (Fig. b) . Radar foreshortening 
occurs in areas with topographic relief when radar waves are reflected 
from sloping surfaces. Radar layover, an extreme case of foreshortening, 
occurs when the top of a feature is imaged before the bottom (Fig. 7). 
Steeply sloping surfaces had to In excluded from the Askja study because 
the r -dar did not portray reflectivity information accurately as a result 
ol geometric distortions in the images. 

Although theories of backscatter ing from terrain are not well 
formulated, it is generally agreed that the terrain parameters which 
affect radar backscatter .ire surface roughness, surface slope, and 
dielectric properties (see the Appendix for a discussion of backscatter 
theory). One of the first empirical geologic investigations using radar 
images was undertaken by Schaber et al. (197b). In their study, radar 


5 




3 


Nt AX IMAG£ fAX 

RANGl RANGl 


Fig. b. Example of radar shadowing (from Sabins, 1978) 



rig. 7. Example of radar foreshortening and layover (from Reeves, 1975) 


Images rf Death Valley were found to show distinctive variations in radar 
backsca:ter that could be correlated with changes in the surface roughness 
of different geologic units. The relationship between radar return power 
and surface roughness for seven backscatter units was evaluated by 
plotting a roughness parameter versus film density across a radar image 
that were converted to relative power. The roughness parameter was the 
sum of the average macrorelief (measured in centimeters), the average 
distance between macrorelief scattering facets (measured in centimeters), 
and the percent of the area covered by micioroughness on the scale of 
A/10, where A is the radar wavelength (Table 2). The significant result 
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Kin. 9. Backseat ter Units 1 — I X determined t'rom radar images of Askja 
Caldera (Unit mt represents pa lagan i tic breccias of 
Dvngjuf joll massif; from Mai In et al., 1978.) 


Nine back:.*, a t ter units were distinguished on the Iceland radar 
images (Fig. 9). Roughness parameters were estimated for six of the 
baekseatter units using data obtained t rom ground photographs (Table 1). 
Tin average vertical macrorel lef , the average distance between macro- 
relief scattering facets, and percent area covered by microroughness were 
estimated from ground photographs which included scales. F.xamples of 
some ol the photos used are shown in Fig. 10. Addltl nal ground photo- 
graphs are presented in Section V. 


Fable ). Roughness parameters estimated 1 rom ground photographs and 
average DN values tor Baekseatter Units 1 —V I in Askja 
region (Units V11-1X were excluded because of insufficient 
data . ) 


h.ck.c.it ter 
unit 

A 

Vert leal 
mac rore 1 let , 

Cl 

B 

Average distance between 
nu it rore l let l acet • , 

cm 

C 

Percent area covered 

: t . 

(A+i*C) 

Roughnc*-* 

parameter 

A*, t.l, . 

DN 

Ke l.tt ftve 
power , 

d» 

I 

15 

10 

40 

115 

190 

-:.i 

11 

3 

10 

35 

40 

m 

-4.0 

111 

2b 

15 

100 

140 

Jlo 

-l.i 

IV 

20 

2b 

40 

135 

2 JO 

0.0 

V 

30 

30 

4 0 

100 

16} 

-3.4 

VI 

b 

b 

5 

15 

no 
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Fig. 10. Examples of ground photos obtained in Askja region (Scale 
marker * 25 cm; photographs by M.C. Mai in.) 

Figures 11 through .’5 are plots of digital number (I)N) versus 
sample location along a given line. The DN values are linearly related 
to power received by the radar since the data were scanner corrected. A 
linear transformation was apnlted to the DN values in order to make them 
proportional to relative powers list'd bv Schaber et al . (Id7f>). A graph 
of roughness parameters versus relative power for the Iceland data is 
shown in Fig. 15. 

The brightness values used in the Death Valley study were determined 
at a constant look angle along one scan line. 'he Iceland scan lines 
were at varying look angles because of the relative locations of the 
backscatter units. The brightness values of two units that span the 
width of the image were compared at varying look angles. No significant 
variation of DN values with look angles could be seen within the accuracy 
of measurements. The lack of statistical data probably results in the 
greatest errors. In the Iceland study, roughness parameters for the 
backscatter units are based on stereograph ic ground photos and field 
observations, while the roughness parameters In Death Valley represent 
averages of many field measurements. 

With these limitations, the graph of roughness parameter versus 
relative power for backscatter units in the Askja region shows trends 
similar to those found in the Death Valley study. A break in slope 
occurs between Unit V (average surface height * id cm) and Unit II 
(average surface height = 3 cm). There are insufficient backscatter ing 
units in the Askja region to determine the breakpoint more precisely. 
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Fig. 12. DN values at sample locations along line 175 
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Fig. 13 . DN values at sample locations along line 300 
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Fin. 15. Graph of roughness parameter versus relative return power for 
Backscatter Units I-VI in the Askja Caldera test site 


Although the Death Valley and Iceland images demonstrate that changes 
in surface roughness can be correlated with variations in grey tones on 
a single polarization image, multiple polarization images must be used 
to determine the overall roughness spectrum of a surface (Daily et al., 
1978 ). 
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SECTION IV 

COMPUTER I MACK PROCESSING 


Digital processing of the Iceland linages was performed using the 
facilities of tiie Image Processing Laboratory (IPl.) at the Jet Propulsion 
Laboratory. The central computer at IP1. Is an IBM 360/65. The software 
system manipulates digitized data and creates a processed image that can 
be converted into a photograph. Computer processing of the Iceland 
Images involved many of the techniques described by Daily et al. (1978), 
Including digitization, scanner correction, registration, and differenc- 
ing and ratioing. 

Radar images were digitized from image film using a high-resolution 
microdensitometer that recorded the amount of light transmitted through 
image film transparencies at selected locations in a rectangular grid. 

In this process, a digital number that corresponds to a grey level in 
the original image is assigned to each picture element. The DN values 
range from 0 (black) to 256 (white). The spacing at which the trans- 
parency is scanned determines the pixel size. For the Iceland data, 
film transparencies were scanned using an aperture that recorded the 
average film density over 40-iim squares spaced 40 ,jm apart, which makes 
each pixel correspond to 20 x 20 m on the ground. 

The power of the radar return is a more linear function of film 
density than of the transmission recorded as DN values. Because the IPL 
film recording hardware exposes film so that film density is linearly 
related to DN, it was necessary to perform an intensity transformation 
on the digitized radar data in order to achieve linearity with the film 
density of the output. 

In order to analyze multiple data sets of the same areas on the 
ground, digital images were spatially transformed to bring corresponding 
surface features to the same location in each image. This registration 
was performed by geometrically transforming one image so that features in 
it overlaid the corresponding features in the reference image. The first 
step in the registration process was to manually select features which 
could be identified in each image (tiepoints). A polynomial that related 
the tiepoint locations on each image was then calculated. The coeffi- 
cients of the polynomial were adjusted to minimize the squared distance 
between the manually selected tiepoints and the locations determined by 
^he polynomial. in many cases, it was difficult to find features that 
could be precisely located in each image. The discrepancies between tie- 
point locations that were measured and those determined by the calculated 
polynomial were checked in order to identify incorrect tiepoints and to 
assess the accuracy with which the computed polynomial modeled the actual 
spatial difference between the two images. In addition, the order of the 
polynomial was varied and discrepancies were reevaluated to determine the 
form of the polynomial that best fit the tiepoints. The polynomial was 
evaluated at control points in a square grid in the reference (master) 
image. Corresponding locations on the other (slave) image were deter- 
mined using bilinear interpolation based upon the four ambient control 
points. Geometrically transformed slave images were generated by assign- 
ing the estimated location in the slave image to each master image 
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location. Figure 4 la a 1 ike- polar ized (IIH) master image. Figure 5 is 
.« geometrically transformed cross-polar ized (HV) image. 

Difference and ratio Images were generated by computing the diff« 
ence or ratio of DN values of corresponding samples of registered imagt 
The differences were displayed with zero difference translated to a DN 
of 128 (intermediate grey). A ratio etjual to one is also displayed as 
Intermediate grey (DN • 50 in this case). Figure lb shows an example of 
an image which is the result of ratiolng a like-polarized (111!) radar 
image to a cross-polarized (HV) radar image. Figure 17 is an example of 
an image which is the result of subtracting an HV image from an HH image. 
The digital form of the ratioed and diffe. fenced data is presented graphi- 
cally in Figures 18 and 11, which are plots of Intensities (DN values) 
along a given line versus sample number (horizontal location). 
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Fig. 16. Polarization ratio image (HH/HV) of Askja radar images 






Fig. 17. Polarization difference picture (HH-HV) o f Askja radar images 
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Fit;. 18. DN values at sample locations along line 525 in polarization 
ratio image (HH/HV) 
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SECTION V 

RADAR RESPONSE TO SURFACE MATERIALS IN THE ASK I A CALDERA REGION 


The overall roughness spectra of the radar backseat ter units in the 
Askja region were predicted from mul t ipy lar ized radar images based on 
studies in Death Valley (see, e.g., Schaber et al. t 1976; Daily et al., 
1978; and the Appendix). A field survey of the caldera was undertaken to 
check the accuracy of the preliminary analysis. Areas visited included 
the caldera floor, the southern rim, and a portion of the southern flunk 
of the caldera. The traverse followed the route outlined in Fig. 20. 

Nine areas with mappable differences were chosen to be field checked 
prior to the traverse. These units are Identified in Figs. 9 and 20. 
However, only eight of the nine areas could be field checked due to time 
limitations. Ground photography of both the overall surface character- 
istics and fine-scale detail were acquired for most areas, including 
oblique and vertical stereo pairs. 

A. RADAR UNIT I (Figs. 21 and 22) 

Unit I is characterized by a high radar reflectivity and a low 
polarization ratio (HH/HV) and difference (HH-VH) . Based on Bragg and 
Rayleigh scattering models (see the Appendix), this unit was Interpreted 
to have a surface of variable roughness with a predominance of scatterers 
greater than 3 cm in diameter. Field observations showed Unit I to con- 
sist of ash and pumice from the 1875 explosive eruption of the crater 
Vlti. In less than 12 hours, the crater, then only 100 m wide, produced 
2 cm“ of rhyolitic pumice. The 1875 eruption was the only eruption that 
produced rhyolitic pumice in the DyngjufjUll area (Wood et al., 1977). 

The surface of Unit 1 is littered with pumice blocks ranging from about 
1 to 40 cm, which is consistent with predictions. 

B. RADAR UNIT II (Figs. 23 and 24) 

Unit II is characterized by an extremely low radar reflectivity. 

This low reflectivity and a polarization ratio of *1 (shown as inter- 
mediate grey on Fig. 6) suggested a smooth surface with a small range 
in surface roughness. Field observations showed that Unit 11 consisted 
of ash and pumice from the 1875 Vi eruption. In some areas, deflation 
and dissection expose extensive deposits of finely bedded ash. The 
surface of these deposits consists of indurated ash less than 0.5 cm and 
numerous clasts of pumice approximately 0.5 to 6 cm in diameter. Defla- 
tion has left cobbles supported on asli pedestals. 

The surface of Unit II was smoother than predicted from reflectivity 
data. This observation and field inspection indicate that moisture may 
be contributing to the radar signature by increasing the radar signal. 

The presence of moisture could also increase the depolarized return more 
than the like-polarized return because of the greater sensitivity of HV 
images to the surface dielectric constant (Daily et al.,1978). Although 
the depolarized return from Unit II is slightly brighter than the like- 
polarized return, this seems to be a minor effect. The etfect ol water 
on the reflectivity cannot he determined accurately since the amount of 
penetration of the radar and the resulting decrease in volume scattering 
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Fi 20. Aerial photograph and map of Radar Backscatter Units 1-IX 
showing traverse (dashed line) and ground photo locations 
(Unit mt represents palagonitic breccias of the 
DyngjufjUll massif; from Malin et al., 1978.) 
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Radar I'n i t II (looking wi st) showing .ish uiul pumice from IH7‘i 
Viti erupt Ion (Pumice clasts range from 0. r j to b cm; 
photograph by M.l . Mai in.) 


tMoseup ol Radar Unit II (photograph by M.C. Matin) 





cannot be ascertained with single-frequency, dual-polarization data. 

C. RADAR UNIT III (Figs. 25 and 2b) 

Unit 111 could be easily distinguished on both the idar image and 
aerial photograph, and it is correlated with a geologic unit called 
Myvetn ingahraun (Bemmelen and Rutten, 1455). l'h’s unit clearly showed 
1 low features. It was interpreted as an aa flow since the radar 
return was strong and the direct to cross-polar izat ion ratio was low. 

My vi tn Ingahraun, which is located in the southwest corner of Askja Caldera 
along the shore of Oskjuvatn, flowed from a north-south fissure in the 
caldera floor near Suderskard in November ’422. The flow has an area of 
1.75 km-’ and an average thickness of 5 m. Field observations showed that 
Unit 111 does have an aa texture (apalhraun, in Icelandic), consisting of 
broken slabs ranging in si/.e from a tew centimeters to a few meters. 

Several areas of underlying Unit 11 ash and pumice show through 
Myvetn inganraun as a result oi incomplete flooding of lava. These win- 
dows are dark patches in the flow on the radar images due to a difference 
in surface roughness, and they are light patches on the aerial photo due 
to a difference in albedo between Units 11 and III. 

D. RADAR UNITS IV AND V (Figs. 27 and 28) 

Unit IV was the second-brightest radar return seen in the Askja area. 
It was brighter than Unit 111, which was a very rough aa flow. The low 
polarization ratio and difierence values also Indicated a greater range 
in roughness than for Unit 111. The reflectivity and polarization ratio 
of Unit V suggested a surface somewhat smoother than those of Units 111 and 
IV. Field inspection revealed both lavas to be the same lava flow. Both 
are lichen-covered aa, but Unit V is partially covered witli volcanic ash. 
The extremely high radar return from Unit IV may result in part from 
the surface morphology ot a hummocky, steep-walled r idge-and-trough 
topography with spacings of 2-3 m and relief of 1 in. In addition, the 
lichen cover may contribute to the greater reflectivity because of its 
higher dielectric constant. The low direct to cross-polarization ratio 
may also be partially caused by the relatively high dielectric constant 
of the lichen enhancing tin HV return. The effect of lichen cover 
cannot be determined accuratelv. Because the dielectric constant 
increases witli an increase in moisture content, the same problem exists 
as in Unit II, that dielectric effects cannot be separated from roughness 
effects on single-frequency, dual-polarization images. 

The return from Unit V is lower than that from Unit IV because of 
the increase in ash cover and the resultant reduction in surface roughness. 
A slight albedo clfference can be seen between Units IV and V on the 
aerial photo that can be correlated with a difference in surface roughness 
on Ihe radar images. 

E. RADAR UNIT VI (Figs. 29 and 30) 

Unit VI is characterized by the lowest reflectivity of all the 
surfaces studied. This unit was interpreted to be smooth. Field 
observations showed that Unit VI consists of a cobble and pebble pave- 
ment cemented by ash, which results in only centimeter-scale relief. 
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Closeup of Rada r Unit 111 (photORraph by M.C. Mai in) 
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Radar Unit IV 
with hummocky 


(looking east) showing lichen-covered 
surface morphology (photograph by M.U 


aa lava 
Mai in) 



Fig. 28. 


K 1 lar Un i t V 
(Same age as 


( looking 
Unit IV; 


west) showing 
photograph by 


ash-cove red aa 
M.C. Mai in.) 


lava 
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Fig. 29. Radar Unit VI (looking southwest) showing cobble and pebble 
pavement cemented by ash (photograph by M.C. Mai in ) 



Fig. 30. Closeup of Radar Unit VI (photograph by M.C. Malin) 
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l'li i * surface produced. a near-specular reflection awa\ t rum the radar. 

Hu' signal returned to the aptcnna Is low in both polarisation modes 
becauHe t he roughness elements are much smaller than t he wavelength of 
the radar. Unit VI can easily be distinguished on the aerial photo In 
its high albedo. 

F. RADAR UN 1 T VII (Figs. 11 and 12) 

The lavas outside the caldera rim to the northeast and east are 
difficult to differentiate on the radar Images. Unit VII has indistinct 
boundaries on the radar images. The overall reflectivity Indicated a 
surface of up to 25 cm relief and a wide range in surface roughness. In 
the tield, the unit consists of aa lava and ash cover. Approximately 20 
of the surface consists of cobble-sized fragments. 

I!. RADAR UNIT VIII (Fig. 33) 

Unit VIII, which occupies the floor of the caldera, is bright in 
both the like- and cioss-polar ized images. The low polarization r.it io 
and low difference value for this unit indicated a wide range in surface 
roughness. In the field, these lavas are often comparable to those ol 
Unit 111 in roughness, although portions have significant ah cove: l rom 
the 1875 eruption. A few outcrops of Unit VI 11 occur along he edge of 
Oskjuvatn. These are completely separated from tin main body of Unit 
VI II by the pumice and ash of Unit 11 that overlies Unit VII 1. 

H. RADAR UNIT IX (no figure available) 

Unit IX displays the best example of lava flow morphology seen in 
the Ask.) a images. The high reflectivity indicated a rough aa surface. 
The alignment of flows is evidence of a linear vent. The highest ref lee 
tivitv in the Askja area and the lowest direct to cross-polarized rat io 
occur at the extreme south end of the flow. This area was not visited 
during the field study because ot logistical problems. Unit IX can be 
distinguished on the aerial photograph and can be correlated with a 
geologic unit. According to Thorarinsson and Sigvaldason (l l »n2), this 
flow erupted between 1^24 and ld2 l >. Aerial photographs and personal 
accounts ot visits to the flow by other investigators suggest a strong, 
similarity between Unit IX and Unit III. 

I. CONCLUSIONS 

In the Askja test site, radar provided additional Information to 
aerial photography since differences in grey tones could be correlated 
with variations in surface roughness. For example, a slight albedo 
difference can be seen between Units IV and V on the aerial photograph 
th.it can be correlated with a difference in surface roughness on the 
radar Images. 


There was good agreement between predicted surface roughness ot 
backscatter units based on computer-enhanced like- and cross-polar ized 
radar images and surface roughness observed in the field. In some cases 
variations in surface roughness could be correlated with previously 
mapped geologic units. The ability to differentiate among lithologies 
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Fig. 31. 


Radar Unit VII 
(Clasts range 


(looking west) showing oa lava and ash cover 
from 10 to 75 cm; photograph by M.C. Malin.) 



Fig. 32. Closeup of Radar unit VII (photograph by M.C. Malin) 
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Rig. 33. Radar Unit VIII overlain by Unit II (looking southwest; 
photograph by M.C. Halin' 



using radar is particularly valuable in areas where Held work is 
impossible, on earth as well as other planets. For example, a good 
understanding of radar baekscatter will aid in interpreting t he 
geology 01 Venus from images returned by the proposed Venus Orbit in, 
Imaging Radar (VOIR) . 
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SECTION VI 
DISCUSSION 


Although single-frequency, dual-polarization images used in the 
Iceland study provided information regarding the average heights of 
surface irregularities and the overall roughness spectra of several 
geologic units, additional frequency and polarization modes would make 
i ossible to discriminate among more roughness units. Multifrequency 
imaging would make it possible to differentiate among more units since 
the transition between diffuse and specular reflection (see the Appendix) 
is wavelength-dependent, according to Rayleigh scattering models. Multi- 
v trequency imaging could also be used to determine the contribution of 

suhsurface multiple scattering because the ratio of a cross-polarized 
image to a direct polarized image is constant as a function of frequency 
when multiple scattering is negligible (Daily et al., 1978). Additional 
polarization modes would make it possible to differentiate among more 
units since simultaneous analysis of HH/VV images and HH, VV, and HV 
images should allow separation of roughness effects from topographic and 
dielectric effects (Daily et al., 1978). The radar interpretation 
techniques, such as those used in the Death Valley and Iceland studies, 
are still being developed. Future studies will concentrate on relating 
surface roughness data to specific rock types. This would involve corre- 
lating roughness spectra derived from mul t ipoiar izat ion and mu It l frequency 
radar data with the texture of primary surfaces as well as texture pro- 
duced by physical weathering of the rock. 
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APPENDIX 


THEORETIC/.: BASIS FOR RACKSCATTKR 1 NG STUDIES 


The average power received by a radar antenna (P^) is given by the 

Radar Equation: „ . 

1 2 2 * 

P-C/\‘o A 

P m — 

K .._v3„4 


(4n) J R 


where 


P ( . * transmitted power 

C “gain of the antenna 

0 

■ normalized backseat ter (fraction of power reradiated from a 
point ) 

A • area of illumination 

R • distance between antenna and scatterer 
X ■ radar wavelength 


The normalized backscatter is a function of radar system properties, 
including polarization, depression angle, and wavelength, and terrain 
parameters, including dielectric properties, surface roughness, and 
terrain slope. The intensity of a radar return is determined by the 
same parameters when all other terms of the radar equation are held 
constant . 


1. Depression Angle 

Depression angle is the angle at which a radar pulse is sent out 
i rom the aircraft, measured from horizontal. It is the complement of 
the angle oi Incidence. Look angle (■ incidence angle) is measured from 
vertical. Therefore, a depression angle of 90 deg and a look angle of 
0 deg are both perpendicular to the surface. The transition between 
■pecular and diffuse reflection from a surface is a function oi depression 
angle and surface roughness. The Rayleigh criterion predicts that 
scatterers with surface heights greater than \/(8 sin 0), where is the 
radar wavelength and ti is the angle of incidence, will result in diffuse 
reflection (Fig. A-l). The magnitude of the normalized backscatter 
coefficient generally increases with decreasing incidence angle. However, 
it becomes independent of incidence angle if the surface is rough with 
respect to radar wavelength (see Fig. A-2). 

2. Wavelength 

Variations in the radar return that are wavelength-dependent are 
directly related to surface roughness and dielectric properties. The 
Rayleigh criterion predicts that a given surface will appear rougher . t 
shorter wavelengths. Attenuation of the radar signal is a function of 
radar wavelength and the complex dielectric constant of a material. The 
depth of penetration of radar below the surface (skin depth) varies 
directly with radar wavelength and inversely witli the imaginary part of 
the dielectric constant. 
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Fig. A-l. Geometry of Rayleigh criterion (h ■ surface height, 
0 ■ angle of incidence; from Srhaber et al., 1978.) 
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lb) ROUGH SURFACE WITH DIFFUSE SCATTERING 
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DEPRESSION ANGLE, d*.j 

(c i RETURN INTENSITY AS A FUNCTION OF DEPRESSION ANGIE 


Fig. A-2. Radar returns from smooth and rough surfaces as a function 
of depression angle (from Sabins, 1978) 
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Polar izat ion 


A like-polarized return is recorded when a horizontally polarized 
radar signal is transmitted and a horizontal 1*. polarized signal la 
received or when a vertically polarized signal is transmitted and a 
vertically polarized signal is received. These modes are called HH and 
W, respectively. A cross-polarized or depolarized return is recorded 
when a horizontally polarized signal is transmitted and a vertically 
polarized image is received or when a vertically polarized signal is 
transmitted and a horizontally polarized image is received. These modes 
are called HV or VH. 

Different types of scatterers will produce different types and 
degrees of depolarization. It is possible to infer what type of scat- 
ters are responsible for an observed type of depolarization (Beckmann, 
1968). Two mechanisms which cause depolarization are (1) multiple 
scattering as a result of target surface roughness and (2) volume scat- 
tering due to nonhomogeneities within the skin-depth of the target sur- 
face (Sabins, 1978). 

4. Terrain Parameters 

The backscatter cross-section of a slightly rough surface has been 
described by Valenzuela (1967 and 1968). The expressions for the dir- 
ect polarization case (HH) and the cross-polarization case (HV) are: 


°HH ’ AU( VV + CI 


0 m> * DW(K ,K ) + F.I 
HV x y 


A A 


2 2 


2 2 
•- r t \ 


A = 4*8 cot 6j(a" cos" 6T^ + sin 6T^)‘ 


C * 2*6^ cot* 6^ sin" 0 sin" 26T" 


D = "8* cot* sin" 26(T., - T^)" 


8 4 •? 7 2 2 2 

E * 2tr8 cot 6^(a" cos" 6 - sin" 6)"T", 
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1'iior*' 

tilt angle perpendlcu lar to the plane ol Incidence* 
0 ■ Look angle 
0. = actual Incidence angle 
a = sin 0 

r " dielectric constant 
Uj = sin 6j 
6 = 2 ti/X 

* = radar wavelength 

w(p,g) ■ energy spectral density of the surface 
For the geometry of the equations, see Fig. A-l. 
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I' i >♦. . A- l. lleomet rv I i'i b.liksrut t rr eiju.lt it'll s 
Uron Daily ft .il., l l »7fO 

Hus. express it’iis i .111 bo iist-il it* lu- Ip uiulf rst .uni tin- behnrloi »* • 
tlu- .1 i t Iti flit 1.. ifksf.it t i'i fi fss-sfi l ifiis .is .1 I line tii'll i'i suit. iii- ill. ii 
art or 1st ies IDullv ft .il., 

i. I In- Kllivt i'i Surture K>' ij’.lnifss 

rin- term W(K N ,K V ) for re f ponds to tin- siin-.lo l»i -»>\>*. sf.it If r inn l » «»m 
tlu- spertr.il f omponent wliifli s.it isl if-, tlu- Hr. inn i-oiiilit ion, \ - ' t •’ 
tM , wlifii- A sui t. iff roughm- ss w.ivi- If until. V radar w.ivol ennt Ii, • • n»l 

look .inn If . On 1 v oin- sportral foiii|'oiu-iit ot tin- surl.iff sat isl’ If s the 
Hr. inn foinlit ioi- tor a nivett wavi- 1 flint Ii anil took .ui>-,lo. lor ox. imp lo, .it 
a Look unnle ft * don. a JS-fm r.nl.u will ho most sons it ivi- to a suit. no 
i oun.hnoss ot approx iin.it f 1 v IS fin. I’livs 1 f.i 1 I v , tin- Hr. inn i' » ' nil it ion rot 
i o .-.points to tin- f.iso win- r i- tlu- projeetod wave vi-ftor is oipi.il to halt 
tlu- wavo voftor ot a rortaln sport ral fomponont ot tlu- surl.iff. I'liis 
if suits in a rosonancf ottfft and a stronn rot urn. 

I In- term 1 is n-latoi! to a double Hi ann rosonanee ottfft, wln-io tlu 
r.ul. 'i wave coml'inat ions eon 1 .1 rosult in ri-sonanff in tills f.iso. Hits 
implies ilia’. I is related to tlu- overall rounluiess spectrum and W is 
related to oniv one eomponent ot that spertrum. lor small tilt angles, 

HU is proport ional to W to tin- l irst ordet ot approxim.it ion. Iliis 
imp lies that imj is also proport ional to tlu- spi-rlial fompoin-ut ot t lie 
surl.iff that sat 1st ies the Hr. inn fond l lion and 'py is proportional to 
tin- ovoi ill rounluiess spertrum. 
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ihe lll.it of Surl.m Slope 


b. 


The effect of surf ace slope is related to both tilt angle perpendic- 
ular to the plane of Incidence and to look angle, since .1 change in look 
angle Is similar to a change in slope In the plane of Incidence. For 
small tilt angles perpendicular to the plane of incidence, A (and up) 
is proportional to cos-’ , while E (and Op) is proportional to cos 4 
Tills impLi* s that py) is more sensitive to than up; Opp becomes sen- 
sitive to surface slope only at 0 larger than about fib deg (see Fig. A-4). 
1 he comparison of the behavior of direct and cross-polarized images as 
a function of look angle cannot be generalized because of the strong 
dependence on W and 1 . 




0. <f|. 


Fig. A- -♦ . behavior of terms A (1111) and F. (HV) in backsactter equation 
(5 » tilt angle perpendicular to incident beam, 

- dielectric constant, 0. “ incidence angle; from 
Daily et al . , 1978) 


< 1 . Hie Effect of Dielectric Constant 

The dielectric constant increases almost linearly with increasing 
moisture content. The backscatter in both 1111 and HV images increases 
with increasing dielectric constant because the reflectivity is higher 
for larger dielectric constant discontinuities at the reflecting surface. 
In the case of large dielectric constant and =0, opp is independent of 
and is proportional to dielectric constant. 
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